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Closed-form expressions are obtained for the impulse response of 
(jraded-index fibers whose relative permittivity is a homogeneous function 
of the two transverse coordinates x, y, and for the impulse width in (jraded- 
index fibers whose profile departs slightly, but otherwise arbitrarily, from 
a square law. The in homogeneous dispersion of the material is taken into 
account. Pulse broadening can be reduced by a factor of 12 from the value 
obtained for square-law fibers. Simple expressions are found for the 
acceptance of highly oversized fibers. 

I. INTRODUCTION 

Light-emitting diodes supply their optical power in a time and space 
incoherent form. The line width is typically of the order of 200 A, 
and the radiation is approximately lambertian with an emissive area 
of the order of 50 X 50 ^m. Time and space incoherent optical pulses 
can be transmitted by oversized optical fibers. However, optical pulses 
propagating in such fibers tend to broaden as they travel. This is in 
part due to the nonzero line width of the source and the dispersion 
(d 2 k/du 2 ) of the fiber material. The other cause of pulse broadening 
is associated with the fact that the time of flight of a pulse along a 
ray depends on the ray trajectory. Pulses traveling along axial rays 
usually go faster than pulses traveling along rays of large amplitude. 
Because both types of rays are excited by spatially incoherent sources, 
the difference in axial group velocity causes a broadening of the input 
pulse. In the main text of this paper, we assume that the carrier is 
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monochromatic and that the spatial distribution of the rays is time- 
invariant. This is the case, for instance, when the source is an injection 
laser that oscillates simultaneously on many transverse modes. The 
difference in frequency between these various transverse modes can 
usually be neglected. 

It was first pointed out by Kompfner 1 that pulse broadening in 
step-index fibers could be drastically reduced by introducing ray 
equalizers at various locations along the fiber. The role of ray equalizers 
is to exchange fast and slow rays. A possible implementation of this 
idea is shown in Fig. 1 together with the calculated impulse response 
for uncorrected and corrected step-index fibers. 2 Because natural 
mode mixing appears to be very small in the most recently made 
optical fibers, ray converters may be practical. They have not been 
experimented with, however, and we shall therefore restrict ourselves 
to uniform, uncorrected fibers. 

Important results concerning the broadening of spatially incoherent 
optical pulses in graded-index fibers have already been reported. In 
Refs. 3 to 9, the difference in group velocity between the various 
modes that can propagate in step-index and graded-index fibers has 
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Fig. 1 — Ray converter that minimizes pulse broadening in step-index fibers, (a) 
Angular spread of a step-index fiber, (b) Optical arrangement with confocal lenses. 
The first and last lenses are unconventional, (c) Calculated impulse response for 
uncorrected [P(0 = 1 within the pulse] and corrected step-index fibers (from Ref. 2). 
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been evaluated. The impulse response is obtained by adding the con- 
tribution of each mode, under the assumption that all modes are 
equally excited by the source. The calculation of the group velocities 
can be simplified with the help of the W.K.B. approximation (see, 
in particular, Ref. 4). 

Let us now describe an alternative ray-optics method. The time of 
flight of a pulse along a ray is first evaluated according to the laws of 
geometrical optics. A ray can be defined by the point x, y where it 
intersects the input plane of the fiber (plane 2 = 0), and by the 
transverse components, k x , k y of the wave vector k. k is, by definition, 
directed along the ray and has magnitude (27r/X )», where X denotes 
the free-space wavelength and n the refractive index of the fiber 
material, usually a function of x and y. Thus, the time of flight of a 
pulse (at a fixed carrier frequency) is, in general, a function of the 
four parameters x, y, k x , k y . These four parameters can be considered 
the components of a four-vector p, in the so-called phase space. The 
impulse response is subsequently obtained by assuming that the 
density of rays is equal to (2tt)~ 2 in the phase space. In other words, we 
assume that the number of rays whose points of intersection with 
the input plane are between x, x + dx and y, y + dy, and whose 
direction is defined by values of k X) k y between k x , k x -+- dk x and k y , 
k„ + dk y , is equal to dxdy dk x dk y /(2ir)' 1 . The total power transmitted 
is the acceptance (or number of modes) of the fiber. This is the power 
transmitted for a source of luminance unity (see, for example, Ref. 10). 

The approach used in Refs. 11 to 13 is based on the conventional 
ray equations. We have shown in Refs. 14 and 15 that it brings a 
considerable simplification to write the ray equations in the Hamil- 
tonian form. The relationship between the ray-optics method and the 
W.K.B. method becomes more obvious with the Hamiltonian form. 
It can be shown that the W.K.B. method and the ray-optics method 
are essentially identical. 14 

An important difference, however, should be noted. In the W.K.B. 
method, modes whose axial wave number k z is less than the free wave 
number k a in the surrounding medium (or cladding) are assumed to 
leak out so rapidly that they can be ignored. On that basis, the accep- 
tance of a step-index round fiber with radius a, for example, is found 
to be 

N 2 = (k 2 - fc 2 s )a 2 /2 =- V 2 /2. 

The radiation loss of leaky modes can be small in the case of highly 
oversized fibers, however, as was pointed out by Snyder. 16 The ray- 
optics condition is distinctly different: Only those rays are ignored 
whose tangential component of the wave vector at the core-cladding 
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interface [_(k 2 z + k 2 ^, where k v denotes the azimuthal wave number] 
is less than the free wave number k B in the surrounding medium. 
According to ray optics, the acceptance of a step-index fiber is N 2 = V 2 
instead of V 2 /2. The influence of the slightly leaky rays on the im- 
pulse response of fibers has not been observed. This is perhaps because 
high-order modes are more sensitive to irregularities than low-order 
modes. Slightly leaky rays may become important when highly over- 
sized fibers of good quality are fabricated. This is even more so for 
graded-index (e.g., near-square-law) fibers, because the field decays 
exponentially beyond the caustic line, which bounds the ray 
trajectories. 

In most previous works, the effect of inhomogeneous dispersion* on 
quasi-monochromatic pulse broadening was neglected. This effect, 
however, was taken into account for square-law and linear-law graded- 
index fibers in Appendix B of Ref. 14, and by Gambling and Matsuhara 9 
for circularly symmetric modes in square-law fibers perturbed by an 
r 4 term. The result for arbitrary small deviations from square-law was 
given by Arnaud in Ref. 15. Olshansky and Keck 9 first pointed out 
that inhomogeneous dispersion is of great practical importance, at 
least for fibers doped with Ti0 2 . Dispersion for the promising Ge02 
doped fibers is not known at the time of this writing. The variation of 
the loss of that material as a function of doping is likewise unknown. 
If we consider further that the sources used in pulse broadening experi- 
ments are not fully characterized in terms of their distribution in 
phase space, it appears that a precise comparison between theory and 
experiment is difficult at the moment. We shall therefore restrict our- 
selves to the theoretical evaluation of pulse broadening. 

II. GENERAL RESULTS 

The derivations of the general results given in this section appear in 
Appendix A. They follow in a straightforward manner from the Hamil- 
ton equations for pulse trajectories in space-time. 

Fibers are most often characterized by a refractive-index profile: 
n(x, y, a)). However, the quantity that enters directly in the expressions 
for pulse broadening is the square of the wave number k 2 (x,y,u>) 
= {2ir/\o) 2 n 2 {x, y, o)), where X denotes the wavelength in free space. 
We shall therefore deal directly with k?(x, y, a>). 

Let x(z), y(z) denote a ray trajectory. Assuming that the fiber is 
time-invariant and uniform and that the material is isotropic, we ob- 



* Inhomogeneous dispersion refers to the spatial variations of the ratio of the local 
phase to group velocities in the material. This parameter should not be confused with 
the parameter d 2 k/du 2 , usually refered to as "material dispersion." The latter is 
important only for broadband carriers. 
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tain from the ray equations the following differential equation (see 
Appendix A) : 

\k*d?{X + Y)/dz i = fc 2 - k* + XdW/dX + Ydk*/dY, (1) 

where we have set, for convenience, X = x 2 (z), Y = y 2 (z). The quantity 
k z in (1) denotes the axial (z) component of the wave vector k and 
is a constant of motion. In other words, k 2 remains the same along 
any given ray. In a wave theory, k z corresponds to the propagation 
constant of a mode (sometimes denoted /3). Note that, in spite of the 
fact that we are using the language of wave optics, the theory given 
in this paper is based strictly on ray optics, except when we impose 
the condition k t > k, to make contact with previous results. 

It follows from the space-time ray equations that the ratio of the 
time of flight of a pulse along a ray to the corresponding time on axis 
is (see Appendix A) 

t = {k<,/k s ){dkydcf)/{dkl/dtf), (2a) 

where A- = k(0, 0, to). The sign ( > denotes an average over a ray 
period. For any function a(x, y, o>), we have denned 

(a(x, y, «)) = Z- J* a\jc(z), y(z), t^dz, (2b) 

where x(z), y(z) denotes a particular ray trajectory and Z the ray 
period. If the ray trajectory is not periodic, (a) should be understood 
as the limit of the right-hand side of (2b) when Z — -> <*> . In the special 
case where the inhomogeneous dispersion of the material can be 
neglected, A: is proportional to w and, consequently, dk 2 /dor = /v 2 /u> 2 . 
Thus, (2a) reduces to 

i = (A- 2 )/(A- A- ; ). (2c) 

Finally, if the source of rays has a distribution /(p) in the phase 
space p = \x, y, k x , k u \, the response of the fiber to an input P'{i) is 
(see Appendix A) 

P(t) - f P'[i - /(p)]/(p)5T(p)(rfp). (3) 

The quantity 7'(p) is the transmission of a ray (usually T < 1), and 
(dp) = dxdy dkjdky. In the special case of a uniform lambertian source 
of luminance unity, we have /(p) = l/(2ir) 2 . For simplicity, we can 
assume that !T(p) is unity when the point x, y is within the core cross 
section and the components k It k y of p are within some area to be 
denned later for specific examples and zero outside that area. All the 
subsequent results follow from (1), (2), and (3). 
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III. IMPULSE RESPONSE WHEN k s (x, y) - k% IS A HOMOGENEOUS 
FUNCTION OF x AND y 

Let the differential equation (1) be averaged over a ray period. 
The left-hand side of (1) vanishes because d(X + Y)/dz assumes the 
same values at the ends of the integration interval. (In this integra- 
tion, k z can be considered a constant.) Thus, we have 

(k 2 - t» + XdW/dX + Ydk 2 /dY) = 0. (4) 

Let us further assume that h(X, Y) = k 2 (X, Y) — kl is a homogeneous 
function of degree k in X = x 2 and Y m y 2 . This means that, for any X, 

h(\X, XF) = X K h(X, Y). (5) 

If we differentiate (5) with respect to X and set X = 1, we obtain 

Xdh/dX + Ydh/dY = K h(X, Y). (6) 

Thus, going back to k 2 (x, y, a>), 

Xdk 2 /dX + Ydk 2 /dY = K (k 2 - kl). (7) 

In that case, a simple and general expression for the relative delay 
in the absence of material dispersion is readily obtained from (2c), 
(4), and (7), 

* = LWfco) + k(a /*.)]/(i + *)• (s) 

The relative delay t is plotted in Fig. 2 as a function of k z /k Q with k 
as a parameter. This result is applicable, for example, to the index 
profile 

k 2 (x, y) = k 2 -a\x\ - (3\y\, (9) 

where a and denote constants. In that example, k = J. Note that 
the fiber described by (9) is not circularly symmetric, even if a = j8. 
Examples of circularly symmetric fibers that satisfy (5) will be given 
in the next section. 

In almost any z-invariant focusing system, any initial distribution 
eventually reaches a steady state. This steady state in general differs 
from the initial distribution. A lambertian distribution / = constant, 
however, remains lambertian because it is a (trivial) solution of the 
Liouville equation (see Appendix A). Note that the distribution / in 
(3) represents a ray (or pulse) density. If the medium introduces a 
nonuniform attenuation on the rays, the power density 7'(p)/(p) in 
phase space needs to be distinguished from the distribution /. 

A fiber is usually surrounded by a homogeneous material, called 
the cladding, with wave number k,. For fibers that are not highly 
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Fig. 2 — Relative time of flight in a fiber where k 2 (x, y) — k\ is a homogeneous 
function of degree k in x and y[k Q = fc(0, 0)]. For most fibers, k//k is close to unity. 



overmoded, the transmission law 
T(x, y, k x , k y ) = 



if k z > k s 
if k z ^ k, 



(10) 



is often acceptable. Equation (10) says that rays whose axial wave 
number is less than the free wave number in the surrounding medium 
are leaking sufficiently rapidly to be ignored. The distribution / of 
the lambertian source is set equal to 1/ (2ir) 2 so that the luminance is 
unity. In that case, the total power transmitted is the acceptance of 
the fiber. The relative time of flight is, within the present assumptions, 
solely a function of k z . The upper and lower bounds on k z are k(x, y) 
and k 8 , respectively. It remains to express the volume element 
dk x dk u dxdy in (3) as a function of dk z , dx, dy. For given x, y, a constant 
value of k z corresponds to a circle of radius squared k 2 (x, y) — k\ in 
the k x , k v space because k 2 x -\- k 2 v = k 2 {x, y) — k'\. Thus, 



dk x dk u = irdk 2 . 



(11) 



Let us evaluate the acceptance of the fiber. The light acceptance of 
any optical system is, as we have seen, the volume in phase space of 
the accepted rays divided by (2t) 2 . It is also equal to the effective 
number of modes that the system can transmit. If we integrate P{t) 
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from t = — oo to t = + oo in (3), the integral over P' in the integrand 
is unity, and we obtain 

.V 2 = (l/4r) If Lk*(x, y) - k\]dxdy, (12) 

where we have u.sed (10) and (11). Thus, 47riV 2 is the volume enclosed 
by the profile : k 2 (x, y) . For a step-index fiber of any shape with cross- 
section area A, for example, we have from (12) 

N 2 = (A/4*-)(A; 2 - A; 2 ). (13) 

This expression should be multiplied by 2 to take into account the 
two states of polarization. 

The pulse transformation in (3) becomes, using (11), 

P(t) = (1/4tt) I dxdy ! ( ' "' P'[t - t(k z )]dk 2 z . (14) 

Let the input pulse P'(t) be a symbolic 8 function (e.g., a rectangular 
pulse of width At and height At~ l in the limit At— »0). The output 
pulse in (14) becomes 

P(t) = (l/^w)\dk 2 /dt\A(k z ), k z > k„ (15) 

where | dk'i/dt | denotes the absolute value of dk\/dt and A (k z ) denotes 
the cross-section area that satisfies k{x, y) > k z . k z can be expressed 
as a function of the delay t by inverting the relation t(k z ) given earlier. 
We obtain, from (8), 

dk' z 2 /dt = 2(1 + K )k'll - (k/K 2 )-], (16) 

where 

K = k z /k = (1 + K)t/2 ± {[(1 + K )t/2J - k}K (17) 

If k > 1, there is only one value of k' z between k' s = k,/k and 1, for 
any k' s . If 

k' 2 < k < 1, (18) 

there are two values of k' z that need be considered. Their contributions 
to P should be added. If k < k' 2 , there is again only one relevant value 
of k' z . 

Let us consider as an example a (noncircularly symmetric) square- 
law medium 

k 2 (x, y) = A:g(l - fl 2 * 2 - Sl 2 y 2 ), (19) 

where ti x , i2 v denote arbitrary constants. 2ir/Q, x and 2-k/SI v , for small 
x, y, are the periods of ray oscillation in the xz and yz planes, respec- 
tively. The area A (k z ) defined earlier is the interior of an ellipse 

A(k z ) = tt(1 - k' 2 )/a x n u . (20) 
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The impulse response is obtained from (15) and (16) with k = 1, 
and (20), 

p(i) = |*WV(°*°»)« *« > k ' (21) 

| , kz = k a , 

where, from (17), k' t = t — (t 2 — 1)*. Because, in most fibers, k z re- 
mains close to ho, the variation of k' z can be neglected, and the pulse 
response is almost rectangular. 

For a step-index fiber, the area A is the area of the core cross sec- 
tion, and t = k /k z . Thus, the impulse response of a step-index fiber 
with cross-section area A is simply 

P(t) = klA/2irt 3 , 1 < t < k /k a . (22) 

Because, in most fibers, t remains close to unity, the pulse response is 
almost rectangular. 2 The pulse width, however, is considerably larger 
than for square-law fibers, as we shall see in more detail later. 

IV. CIRCULARLY SYMMETRIC FIBERS WITH Jr s - *?, A POWER 
OF THE RADIUS 

Let the wave-number profile be of the form 

k 2 (R,u) =k 2 (co) -kl(u>)R<, (23) 

where R = X -\- Y = r 2 denotes the square of the radius. The relative 
time of flight is, substituting (23) in (2a), 

t = (h/k,)(dk 2 /da>*)/(dkl/du 2 ) 
= (k /k t )(l - e K D K (R<)), (24) 

where we have defined 

e« - kl/k 2 (25) 

D K = kKdkl/d^/kKdkl/cM*). (26) 

D K is a dispersion factor equal to unity in the absence of dispersion. 
Thus, we need to evaluate (R K ). It is interesting that we can do that 
without solving the ray equations. The quantity (R K ) is, of course, 
independent of dispersion, so we may omit the w arguments. 
For circularly symmetric fibers, (1) can be written 

±kid 2 R/dz 2 = d(k 2 R)/dR - ki (27) 

Averaging (27) over a ray period, we obtain 

ki - (d(k 2 R)/dR). (28) 

We have also, directly from (23), 

(A: 2 ) = kl - ki(R') (29a) 
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and, from (28) and (23), 

(k 2 ) = (kl + Kkl)/(1 + k). (29b) 

Equating the two expressions (29a) and (29b) for (k 2 ), we obtain 

u(R') = (l - K 2 )/(l + k), (30) 

where k' z = k z /k . Thus, substituting (R K ) from (30) into (24), the 
relative time of flight is 

t = k'r 1 - DM" 1 - K)/(i + k). (3i) 

In applications, we need k' z as a function of t. Solving (31) for k' z and 
setting D' K = D K /(l + *), we obtain 

k' z = (t/2D' K ) ± l(t/2D' K y + 1 - D;- 1 ]!. (32) 

By differentiating (32), we further obtain 

dk' a a /dt = ZkXpi - (1 - D' K )/k' t 2 ']-\ (33) 

To obtain explicitly the impulse response in (15), we need the area 
A(k z ) defined by k z < k(R). For k(R) in (23), this area is 

A{k z ) = wR(k z ) = x[(l - k' z 2 )/e K J". (34) 

If e, were kept a constant as the parameter k varies, the core radius a, 
defined by k(a) = k 8) would vary. Thus, it is preferable to express e K 
as a function of the core radius a. We have 

e >A = (1 - k'W'/a*, (35) 

where k' s = k s /k . The impulse response is finally obtained from (15), 
(33), (34), and (35); 

P(t) = (*8a»/2)*;[(l - k z 2 )/(l - W/W - (1 - DZK-*l (36) 

The possibly doubled value k' z is expressed as a function of t by (32). 
Thus, a closed-form expression has been obtained for the impulse 
response of a fiber with k 2 — kl a power of r, that takes inhomogeneous 
dispersion into account. 

In the absence of dispersion, we have D' K = 1/(1 + k), and (36) 
reduces to 

P(t) = (*ga»/2)*I(l - AVVU " *W'<1 + «)/(! - Kk z ~ 2 ). (37) 

As indicated in the previous section, there are in general two values 
of k' z that contribute to P. Note that the shape of the impulse response 
does not depend on the core radius a. 

The impulse response P(t) in (37) is shown in Figs. 3 and 4 for 
various values of the parameter k. These curves are essentially the 
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l 

Fig. 3 — Impulse response of a circularly symmetric fiber with k 2 (r) = k% — k\r u 
for a lambertian source and various values of k. The optimum impulse response is for 
k w k'. = 0.9. 




Fig. 4 — Continuation of Fig. 3 for larger values of k. 
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same as those shown in Ref. 4. Figures 3 and 4, however, are much 
more detailed. We have assumed that k,/k = 0.9, that is, An/n = 10%, 
a rather large value. For k = 1 (square-law fiber), the pulse width t 
is 0.0054. For example, if n = 1.45 and the fiber length is 1 km, the 
pulse width is 26 ns. For k = 0.9, however, the corresponding pulse 
width is only 7 ns. We find, in agreement with Ref. 4, that the mini- 
mum pulse width occurs when k = k' s . For a step-index fiber (k — »°°), 
the pulse width would be as large as 630 ns. Note the following detailed 
features on the curves in Figs. 3 and 4. For (0.9) 2 < k < 1, the re- 
sponse starts from infinity because of the minimum in the t(k z ) curve. 
For k = 0.85, P drops suddenly for t ^ 0.998. This is because, at that 
time, the smaller of the two k' z becomes less than 0.9, and is rejected. 
For k = 0.95, the response crosses the t = 1 axis. 

Figure 4 is applicable to larger values of k. We note that, for a 
very large k (step-index fiber), the response is almost rectangular. The 
slow decay in power shown in Fig. 4 would be almost negligible for 
small An/n. 

The effect of inhomogeneous dispersion is shown in Fig. 5. The 
parameter k is kept equal to 0.9 (this is the optimum value in the 
absence of inhomogeneous dispersion), but D K is made to vary in the 
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Fig. 5 — Impulse response for a fiber with k 2 (r) = k% + k\(r 2 )°f for various values 
of the parameter D that expresses the inhomogeneous dispersion of the material. 
D = 1 corresponds to the absence of dispersion. D ^ 1 merely introduces a shift 
in the optimum value of «. 
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neighborhood of unity. These curves have a striking resemblance to 
those in Fig. 3. This means that the effect of inhomogeneous dispersion 
merely consists in shifting the optimum value of k. The impulse re- 
sponse remains essentially the same, at least for K & 1. 

The total pulse power is the acceptance of the fiber, a function of k. 
The acceptance is, in the present case, 

N 2 = ( +CO P(t)dt = (J) J* [k*(R) - tffidR 

= (i)C(fcg - fc 2 .)o 2 - WWV(« + i)] 

= [k/4(k + l)](A-5 - fc 2 ,)a 8 . (38) 

The coefficient in the last expression in (38) is f for step-index fibers 
(k — >oo) and | for square-law fibers. The acceptance given in (38) 
should be multiplied by 2 to account for the two states of polarization. 
The same rule applies to all the expressions given in this paper. It is 
more difficult to obtain the ray-optics acceptance of fibers. The result 
is derived in Appendix B. 

In the next section, we consider fibers whose profile departs slightly, 
but otherwise arbitrarily, from a square law. 

V. NEAR-SQUARE-LAW FIBERS 

Let us rewrite the differential equation (1) for circularly symmetric 

fibers 

\k\d l R/dz* = d(k-R)/dR - k%. (39) 

For square-law fibers 

k-(R) = A-5 - k\R, (40) 

the solution of (39) is 

R{z) = R + (fiij - li/kiy> cos mz/k' z ), (41) 

where 

R Q = ±(k 2 - m/k\ = *(1 - k?)/W (42) 

and ft = ki/k . We have introduced in (41) the axial component of the 
angular momentum (or Bouguer invariant) 

U = xk y — yk x , (43) 

which is the second constant of motion. Let us set 

9 = {U^RmY, (44) 

where Rm denotes the maximum radius squared. Note that, for merid- 
ional rays, = and, for helical rays, 0=1. Equation (41) can be 
written in the convenient form 

R = $R M (l + 0) + $R*i(l - 0) cos (2Qz/K). (45) 
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For later use let us evaluate (R n ), the average of R" over a ray 
period. Using the binomial expansion and the result 

(cos" 1 ) = m!2-"»[(m/2)!]- 2 (46) 

for m even and for m odd, we obtain 

{B ) - n!2 ^/ m _L 2M(n _ m)![(m/2)!]2 - (47) 
In particular, 

(R 2 ) = Rl[(S6 2 + 2d + 3)/8 (48a) 

(R*) = fl2/(l + 0) (50 s - 2d + 5)/16 (48b) 

<ft 4 > = Rl,(35e i + 2O0 3 + 180 2 + 205 + 35)/12S. (48c) 

Let us now show that a closed-form expression can be obtained for 
the times of flight in fibers whose permittivity profiles depart slightly 
from a square law. Inhomogeneous dispersion is taken into account. 
Let the profile be of the form 

k 2 (R) =kl-k\R+ f k 2 n R\ (49) 

We assume that e H R"~ l , n ^ 2, is of the order e « 1, where e„ = kl/kl. 
Substituting (49) in (2a), we obtain (with ft 2 = ei = k\/kl) 

t - £-»(l - Di&ifl) + £ e n D n (R«)), (50) 

H=2 

where we have defined inhomogeneous dispersion factors 

D n = {kldk 2 n /do?)/(kldkl/drf). (51) 

The D n are unity in the absence of inhomogeneous dispersion. Because 
the perturbation is small, (R") in the sum (50) can be replaced by 
the expression (47) applicable to square-law fibers. This approximation 
is not permissible, however, for the term (R) in (50) because this term 
is not small. We need an exact expression for (R). We proceed as in 
the previous section. We first observe that, for k 2 in (49), 

d(k 2 R)/dR = 2A; 2 - k 2 + £ (n - l)fc*fl». (52) 

B— J 

Integrating (39) over a ray period, the left-hand side vanishes and, 
using (52), we obtain an expression for (A; 2 ) that does not involve 
(R) 

(k 2 ) = l(jbi + kl) + * £ (1 - n)k*(R"). (53) 

n=2 

1192 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1975 



We also have, directly from (44), 

<fc 2 > = k 2 - k\(R) + £ k 2 n (R"). (54) 

n-2 

Thus, by comparing (53) and (54), 

hl(R) = Htt ~ kl) + h Z (« + l)*2<fl">. (55) 



Substituting this expression for (R) in (50), we obtain the relative 
time of flight for circularly symmetric near-square-law fibers 

t = fcrMi - Mi - ^ 2 )£i + £ U>* - |(n + i)/)i]€.<fl->|. (56) 

71=2 

Alternatively, J can be expressed in terms of the azimuthal and radial 
mode numbers. The result is given in Appendix C. 

In the absence of inhomogeneous dispersion, (56) reduces to 



.v 
I 

n=2 



t = &T l K + 1 + £ (1 - n)e n (R*)J (57) 



Limiting ourselves to an r* correction to the square-law profile, 
C3 = C4 = . . . =0, and setting e 2 = e, (57) becomes, using (48), 

t = i[l - PM (1 + 0)]-*[2 - p,u(l + 6) - € PM (3d* + 20 + 3)/8] 
W 1 + pS,[(2 - 3e) + (4 - 2e)0 + (2 - 3 C )0 2 ]/16 + 0(p 3 M ) 
p,r = &R„. (58) 

The first two terms in (5S) give sufficient accuracy when pm < 0.01, 
that is, when the total relative change in refractive index An/n ^ p a /2 
is less than 0.005 (p« = fi 2 a 2 ). 

The total pulse width r is the maximum variation of t for < 6 < 1 
and < p < p n . For the square-law fiber [e = in (58)], we obtain 

t = p'i/2 (ray optics). (59) 

It should be noted that, in defining t in (59), we have specified that 
the maximum radius of the ray be less than a for any 8. This condition 
is different from the condition used earlier that k z be larger than k„. 
The ray-optics condition pit < Pa is applicable to highly oversized 
fibers. 

If we now consider the expression in (58) with a correction term in 
r 4 , we find that t — 1 f or meridional rays (0 = 0) when e 2 = f in 
agreement with Ref. 17, where it is shown that all the rays have 
exactly the same optical length when fc 2 (x) = [cosh (a:)] -2 £d 1 — x 2 
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+ (!)* 4 + ■ • •• We also find that t = 1 f or helical rays (0 = 1) when 
62 = 1, in agreement with Ref. 3, where it is shown that helical rays 
have exactly the same optical length when k 2 (r) = (1 + r 2 ) -1 ^ 1 
_ r 2 _|_ ? .4 _|_ ... gy considering all rays whose maximum radius is 
less than a, we find that the minimum t is obtained for e 2 = 0.91. In 
that case, r = 0.046pa- The improvement compared with square-law 
media is therefore as large as 11. If we had imposed instead the wave- 
optics condition k z > k„ the vertical scale in Fig. 6 would be divided 
by (1 + 0) 2 . For e 2 = 0, for example, the wave-optics pulse width is 
Pa/8 instead of p 2 /2 as in (59). With the wave-optics limit, the opti- 
mum value of e 2 turns out to be f instead of 0.91. The improvement 
over the square-law case is only 4, instead of 11. 

Let us now consider the effect of r G terms. Figure 7 shows the varia- 
tion of the pulse width t, defined as the maximum variation of t for 
any < 6 < 1 and any < pjg < 0.002, as a function of e 2 for various 
values of e 3 . The effect of e 3 is essentially to shift the optimum value 
of €2 to lower values. The reduction in pulse width is rather modest. 
Nevertheless, a small improvement is obtained, compared to the case 
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Fig. 6 — Variation of the normalized time of flight for a fiber with k 2 (r) = kl — k]r 2 
+ ej(A:}/A;g)r* in the absence of material dispersion for various values of the parameter 
«2. = corresponds to meridional rays and 6 = 1 to helical rays. «2 has been 
redefined to be dimensionless. 
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Fig. 7 — Variation of the width of the impulse response with « for various values 
of e, for a fiber with k*(r) = k* - ky + ei(k*/kl)r* + e 3 {k\/kQr*. 



where 63 = when 



/b - 2 /fc 2 (r) = 1 - p + 0.615 P 2 + 70p 3 , p a = 0.002 
p = tPf*. 



(60) 



We give only the result when the departures from a square-law 
profile are not circularly symmetric. The free wave number in the fiber 
is now in the form 

k*(x, y, «) = A: 2 (c) - ki^R + £ t k 2 nl (u>)X l Y»- 1 , (61) 

where, as before, X = x 2 , F = y\ R = x* + ?/ 2 = r 2 . The ratio, t, of 
the time of flight along a ray to the corresponding time on-axis is found 
to be 

t - k'-Mi - Id - W + £ £ P>.i - *(" + ^d 

71=1 f=0 

X€ni(X l Y*-*)) t (62) 

where e n( = k 2 nl /k 2 Q and D nJ is defined as D n in (51) with A; n replaced 
by Ar„i. Let us assume that it is permissible to use the sinusoidal rays 
of the square-law medium to evaluate the quantity (X'Y"-'). Because 
the average over one cycle of the product of powers of sinusoidal func- 
tions is known, the relative delay t can be written in closed form. 
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Let the ray trajectory be written 

x(z) — x Q cos (az + <f> x ) (63a) 

y(z) = y cos (az + <t> v ). (63b) 

The coefficient a does not enter in the final result and is henceforth 
omitted. We evaluate 

(X l Y n ~ l ) m (a*y<— ») 

= 4^ (n -"<[cos (z + a )J'[cos (0 + 0„)] 2 <"-'>)- (64) 

It can be shown that 18 

([cos (z + 0*)] 2 Tcos ( 2 + *„)]'<»-'>) 

=H 2 I' (» - r- J d - .) cos *•<* - *» )] 

where 

(ft) S (a-l)!6!' (66) 

Thus, given a ray trajectory, defined by the parameters x , yo, <f> x , 
and 4> y (or, equivalently, by the values of x, y, k x , and k v at the input 
of the fiber), we can evaluate in closed form the quantity (X l Y n ~ l ) 
that enters in formula (62) for the relative time of flight, from (64) to 
(66). 

The above calculation is incomplete for the following reasons. When 
the power law n?(r) of a fiber departs from the exact square law, pro- 
jected ray trajectories in the (xy) transverse plane are precessing 
ellipses.* That is, the principal axes of the near-elliptical trajectories 
slowly rotate as a function of z. This precession is unimportant for 
circularly symmetric fibers. For noncircularly symmetric fibers, how- 
ever, the ellipse precession introduces an averaging effect. Further- 
more, the noncircularly symmetric components of n 2 (r, <p) change the 
eccentricity of the precessing ellipse. The axial component k z of the 
wave vector remains a constant, but the axial component l z of the 
angular momentum varies. Finally, in real fibers, slow (adiabatic) 
changes of the refractive index law along the fiber axis are likely to 
occur that must be taken into account. The twists of the fiber axis 
must also be taken into account. Thus, a realistic assessment of the 
effect of small noncircularly symmetric departures of the index law 



* It is well known in mechanics that the only r u potentials (potential U ~ n 2 ) 
that give closed trajectories are the harmonic potential U(r) ~ n 2 (r) = 1 — r 2 and 
the Newton potential U(r) ~ n % (r) = 1/r. 
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from square law on pulse broadening requires a deeper and more 
intricate analysis than the one given in the present section. However, 
the result in (62) and (65) can be used as a basis for more complete 
analyses. 

VI. CONCLUSION 

From a rather straightforward application of the Hamilton ray 
equations, we have obtained closed-form expressions for the pulse 
width in graded-index fibers when k 2 (x, y) — k% is a homogeneous 
function of x and y, and for fibers whose profile departs slightly, but 
otherwise arbitrarily, from a square law. Inhomogeneous dispersion 
was taken into account. The expressions obtained are exact. The small 
angle (or weakly guiding) approximation need not be made. We have 
also given simple expressions for the wave optics and ray optics 
acceptance of weakly guiding graded-index fibers. 

The algebraic results given should prove more accurate and require 
much less computer time than the straightforward numerical integra- 
tion of time along ray trajectories. We have carried the perturbation 
only to first order in the small parameter e. To obtain more accurate 
results, up to order e 2 , we need a more accurate expression of the ray 
trajectory. This expression can be obtained, for example, by the 
method of strained coordinates. 19 These more accurate expressions are 
probably not needed, however, in most practical cases. 
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APPENDIX A 

The Hamilton Equations of Ray Optics 

The Hamiltonian form of the ray equations are well known in 
mechanics and wave dynamics, 20 and they have also been used fre- 
quently in optics (e.g., Refs. 14, 15, 17, 21, and 22). However, their 
simplicity and power is not always appreciated. The physical difficulty 
is that it is not always recognized that ray momenta and wave vectors 
(or photon momenta) are identical concepts. On the other hand, ray 
momenta (proportional to the wave vectors) need be carefully dis- 
tinguished from mass-carrying momenta (proportional to the group 
velocities). 23 On the mathematical side, we need distinguish a function 
such as k z (x, y) and the value k z assumed by that function. We must 
also be aware that da/dz denotes a total derivative, that is, in the 
present context, the variation of the quantity a along some given ray. 
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If a is a known function of x and y, and x = x(z), y = y(z) denote a 
known ray trajectory, then da/dz = (da/dx)(dx/dz) + (da/dy)(dy/dz) 
can be evaluated explicitly as a function of z. Here again, an arbitrary 
point in space x, y should not be confused with a specific ray trajec- 
tory x = x(z), y = y(z). Unfortunately, it is not possible to go into 
more details here. An excellent reference is Lighthill's paper. 20 A com- 
parison between the W.K.B. method and the Hamilton equations is 
given in Ref. 14. 

Let X = (x, y, z, id) denote a point in space-time (t is time) and 
K = (k x , k ti , k z , ico/c) denote the four-wave vector, with w the angular 
frequency. An arbitrary medium is characterized by a function of K 
and X that we denote 

H(K, X) = 0. (67) 

The Hamilton equations for light pulses X(o-), K(o-) are 

dX/da = dH/dK (68a) 

d,K/d<r = -dH/dX, (68b) 

where a denotes an arbitrary parameter. 

Equations (68a) and (68b) can be considered the basic postulates 
of geometrical optics. From a wave-optics point of view, (68a) follows 
from the requirement that the wave lengths and periods of the waves 
that constitute a wave packet be the same in the direction of a ray. 
Equation (68b) follows from (67), (68a), and the fact that K is the 
gradient of an eikonal function. Thus, in wave optics, the Hamilton 
equations (68) are derived from first principles and need not be 
postulated. 

Let £ denote a point in phase space (k x , k y , cc, x, y, t) at the input 
plane, and £' a point in phase space at the output plane. The optical 
system maps the input phase space into the output phase space, that is, 

* = mr (69) 

It follows from (67) and (68) that the Jacobian of the transformation 
(69) is unity, a result often used in photometry. Equivalently, we can 
say that the determinant of paraxial ray matrices is unity or that the 
ray density in phase space is a constant of motion (Liouville theorem). 
These three statements are obviously equivalent, provided the rays 
are not reflected. 

A source of light that is time and space incoherent is described by 
a distribution S(Q in phase space. Each small volume in phase space 
can be pictured as an optical pulse, provided the dimensions of the 
volume are larger than unity. More precisely, this picture requires 
that AuAt » 1, Ak x Ax^>l, and Ak y Ay^>>l. The detailed structure 
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lA* 



of the pulse is ignored in ray optics. Only the motion of the center is 
considered. 

The transmission 7\ of an optical pulse through the optical system 
is a presumably known function of £ that we denote as 

T x = 7\(0- (70) 

For lossy optical systems, T\ < 1. Because the Jacobian of the trans- 
formation £ — * fj is unity, the output distribution is simply 

«'(£') - S(©Ti(0. (71) 

The power emitted by the source and the power that can be collected 
at the output of the optical system are obtained by integrating S 
(or S') over all variables, except t (or t'). Thus, 

Pit) = J 5(0 (dC) (72a) 

P'(t') = J S'(Z')W), (72b) 

where < = (k x , k y , w, x, y) and <' is similarly denned. The terms (d<) 
and (d{') denote elementary volumes in { and (' spaces, respectively. 
The response of the detector could be described by a function D(£'). 

For simplicity, we do not take the detector response into considera- 
tion. All subsequent results follow in a rather straightforward manner 
from the above results, through a succession of approximations. 

Let us assume that the properties of the fiber do not vary with time. 
This means that the Hamiltonian in (67), the transmission T lt the 
mapping £ — ► £', and the pulse delay do not depend on time. In 
particular, 

t' = t + h(0. (73) 

Sources that are ^-separable, on the other hand, have the property that 

SU) = P(t)F(0- (74) 

That is, the distribution in (-space does not vary with time. For a hot 
tungsten wire whose temperature varies with time, the spatial phase- 
space distribution is almost lambertian at all times, but the frequency 
spectrum (approximately given by the Plank law of radiation) varies 
with time. Thus, (74) is not applicable to that source. For consistency 
with (72), we assume that F(Q is normalized to unity. 
For most sources, we can further assume that 

F(C) = B(«)/(p), (75) 
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where p = (k x , k v , x, y) denotes a point in spatial phase-space. That 
is, we assume that the spatial distribution does not depend on what 
part of the frequency spectrum we are considering. Both ft and / are 
assumed normalized to unity. This ensures that F is normalized to 
unity. When the spectral width of the source is small (e.g., less than 
1 percent, as is the case for light-emitting diodes) and the fiber material 
absorption does not exhibit sharp resonances in that band, we can 
assume that 

Ti(0 = r„(«)r(p) (76) 

and 

*i(0 = *o(«) + t(p). (77) 

For definiteness, we assume that the maximum value of Tq(w) is 
unity, and we define t (co) as the delay experienced by axial pulses. 
We evaluate in the main text t(p)/to at a fixed angular frequency. 
The pulse response is obtained from (71) to (77), 

P'(0 - j Pit' - t (u>) - <(p)]S2(a J )!To((o)/(p)7 1 (p)(dp)rfa, 

= f P"lt' - « (co)]S2( W )7 , o(co)d W , (78) 

where 

p"(n = / Pit" - *(p)]/(p)n P )(d P ). (79) 

In writing (78) we have used the fact that the Jacobian of the trans- 
formation £ — > £' is unity and that dw = dw'. The pulse response is 
the convolution of the pulse response in (79), which we may call the 
quasi-monochromatic pulse response, and the spectral width of the 
source. In most cases, T (co) is a constant. For injection lasers, the quasi- 
monochromatic pulse response is the most important contribution. 

In what follows, we assume that the fiber is uniform and long com- 
pared with the period of ray oscillation and therefore approximately 
z-invariant. Let the Hamiltonian in (67) be written 



H = k x — k z (k x , k y , 03, x, y) = 0. 


(80) 


The Hamilton equations (68) are 




dx/dz = —dk z /dk x 


(81a) 


dy/dz = —dk z /dk v 


(81b) 


dt/dz = dkz/du 


(81c) 


dkjdz = dk z /dx 


(81d) 


dk v /dz = dk z /dy. 


(81e) 
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Let us assume further that the medium is isotropic, that is, 

k\ = k 2 (o>, x, y) -k 2 z - k\. (82) 

Thus, (81a) to (81e) are 

dx/dz = k z /k z (83a) 

dy/dz = kjk z (83b) 

dt/dz = (dk 2 /da)/2k z (83c) 

dkjdz = (dk 2 /dx)/2k z (83d) 

dky/dz = (dk 2 /dy)/2k t . (83e) 

According to (83c), the time of flight of a pulse along a ray for a period 
(period = Z) is obtained by integrating (dk 2 /dco)/2k z from z = to 
z — Z. If fco(w) = k(u, 0, 0) denotes the wave number on axis, the 
time of flight of a pulse along the z axis is similarly obtained by inte- 
grating (dko/dco)/2k . Thus, the ratio of the time of flight of a pulse 
along a ray to the corresponding time on axis is 

t = (k,/k z ){dkyd<j)/{dki/d<f), (84) 

where ( ) denotes an average over a ray period. If the trajectory is not 
periodic, ( ) is understood as a limit for z — > « . When k is proportional 
to a), (84) reduces to 

t = (k 2 )/k k z . (85) 

Let us now observe that, from (83a), (83b), (83d), and (83e), 

%k\d*{X + Y)/dz 2 = k 2 - k\ + XdkydX + Ydk 2 /dY, (86) 

where X = x 2 , Y = y 2 . This is easily verified by carrying out the differ- 
entiations. Equations (86), (84), and (79) (with a slightly different 
notation) are those used in the main text. 

APPENDIX B 

Acceptance of Highly Oversized Fibers 

The acceptance, or effective number of modes transmitted by the 
optical system, is the volume of the accepted rays in phase space 
divided by (2ir) 2 . We have said earlier that, if the fiber is very long, 
all leaky rays are eliminated and the acceptance is simply the volume 
enclosed by the profile k 2 (x, y) divided by 4tt. If the fiber is highly 
oversized, however, many leaky rays {k z < k a ) are not significantly 
attenuated. 16 We need then consider the ray-optics condition that the 
tangential (rather than axial) component of the wave vector be larger 
than k, at the core-cladding interface. The ray-optics acceptance is 
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now evaluated for circularly symmetric fibers. We specify that 

k\ + k% > fci a t r = a, (87a) 

where k v denotes the azimuthal wave number at the interface. We 
also have the condition 

k\ > 0, (87b) 

which is not implied by (87a). In this appendix, we restrict ourselves 
to small differences in refractive index, in which case condition (87b) 
can be ignored. Because of the conservation of h (the axial component 
of the angular momentum), we have 

rk u = ak v (88) 

for a ray with x = r, y = 0, k x , k y , at the input plane, that can reach 
the interface r = a. Thus, condition (87a) is 

fc 2( r) _ k 2 _ k 2 + (f/ a *)kl > kl (89) 

Equation (89) defines an area in the k x , k y plane bounded by an ellipse. 
We have to make sure, however, that rays outside that area do in 
fact reach the interface. This is not necessarily the case. The maxi- 
mum ray radius tm is defined implicitly by 

k 2 x + (1 - ryrfok* = k 2 {r) - W(r M ), (90) 

where r\i is the largest real number that satisfies (90). (The initial 
radius r is considered a constant in the present discussion.) Equation 
(90) shows that the k x , k y that correspond to Tm are contained in an 
ellipse with semi-axes squared k% = k 2 (r) — k 2 (r M ) and k^ = [k 2 (r) 
—k 2 {r M )~\/{\ — ? ,2 /rS/), respectively. If k 2 (r) is never increasing, we 
are sure that k xQ keeps increasing as Tm increases from r to a. We do 
not have any such assurance for k y0 , however. When Tm reaches a, 
there may be acceptable values of k x , k y that are located outside the 
ellipse defined above. For each profile, we need therefore verify that 
klo( r M) never exceeds k 2 yQ {a). We easily verify that this is the case for 
square-law fibers, because 

k 2 y0 = mis - r*)/(l - r 2 /r 2 M ) = kY M (91) 

increases with tm for any r. 

Thus, for square-law fibers at least, we can proceed with the calcu- 
lation of the area of the ellipse defined by (89). This area is 

7r[fc 2 0-) - fc^l - f-/a 2 )-K (92) 

Substituting this result in the general expression for the acceptance 
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factor, we obtain 

N 2 = (1) f ' \V(r) - fc 2 ](l - fVa*)-»dr». (93) 

This expression simplifies if we introduce the variable u = (1 — r'/o 2 )*. 
Equation (93) becomes 

AT* = (a 2 /2) f [fc 2 (u) - fcQdu. (94) 

Jo 

Thus, the ray-optics acceptance of most circularly symmetric fibers is 
half the area enclosed by the curve k 2 (u)a 2 . For a step-index fiber, we 
obtain from (94) 

N 2 = (fcg - h*)a 2 /2 (step-index, ray optics). (95) 

This is twice the wave-optics acceptance. Thus, for step-index fibers, 
the slightly leaky rays carry half the power. Our result agrees with 
that in Ref. 16 for weakly guiding fibers. For a square-law fiber, with 
k(a) = k t , we obtain 

N 2 = (fcg - /c 2 )a 2 /6 (square-law, ray optics). (96) 

In square-law fibers, 25 percent of the total power is carried by slightly 
leaky rays. * 

APPENDIX C 

Impulse response width of near-square law fibers 

When the source distribution is lambertian, all propagating modes 
are equally excited. It is convenient in that case to express the relative 
time of flight t for near-square-law fibers given in (56) as a function of 
the mode numbers (azimuthal number: n = ■ • ■ — 2, —1, 0, 1, 2 ••• 
and radial number: a = 0, 1, 2 • • •) rather than k z and l z . This can be 
done by quantizing the ray trajectories. [If the W.K.B. method is 
used, it is essential to first remove the singularity of the Helmholtz 
equation at r = 0. This is achieved by changing the independent 
variable from r to log (r).] One easily finds that the axial component 
of the ray angular momentum l z is equal to n. Furthermore, we can 
use for k z the well-known expression applicable to square-law media. 
The result (56) is written below as a function of a, n, for the reader's 
convenience. We have 



t(a, M ) = (1 - B)~i (l - hBD, + t/yNy) , 



(97) 



* This is in agreement with a recent result by D. N. Payne. 
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where 

B ^ 2gK\ /2 /K 

g s 2a + | A* | + 1; ImI— abs. val. (/i) 
F y m 7 !2->[D 7 - |(7 + l)D{\KJ '{K*Rl'*) 

N, m (2g)y ± {2-( 7 - m)![(m/2)!]'}-i[l - (v/gW*. (98) 

m-0,2 

The parameters i? 7 = ^, 7 = 0, 1 • • • and Z) 7 , 7 = 1, 2 • • • are ob- 
tained from the square of the wave number: K(R) = k 2 (R) = (co/c) 2 
?i 2 (.ff) of the fiber as a function of R = r 2 , measured at the nominal 
wavelength \ and at a slightly different wavelength, Xq, expanded in 
power series of R as follows 

K(R) = K - K X R + K 2 R 2 + • • ■ (X ) (99) 

K'(R) - Ko - K[R + K' 2 R* +■■■ (X„). 

The D y are obtained from (99) 

D y = K (K' y - Ky)/K y (Ko - K ) . (100) 

If we can neglect the power in the leaky modes, the mode numbers 
a, /1 are restricted by the condition k z > k„, that is, 

B < 1 - K./Ko « 2An/n, (101) 

where K, = k 2 is the square of the cladding wave number. The root- 
mean-square impulse response width is defined as 

a = 5,000[<£ 2 ) - <*) 2 ]* ns/km, (102) 

where ( ) denotes an average taken over all the modes permitted by 
(101). Thus, it is a straightforward matter to evaluate from our ex- 
pression in (56) the root-mean-square width of the impulse response of 
any circularly symmetric near-square law fiber, provided the wave- 
number profile can be measured with sufficient accuracy at two 
wavelengths. 
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